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Current observational bounds on dark energy depend on our assumptions about the curvature of 
the universe. We present a simple and efficient method for incorporating constraints from Cosmic 
Microwave Background (CMB) anisotropy data, and use it to derive constraints on cosmic curvature 
and dark energy density as a free function of cosmic time using current CMB, Type la supernova 
(SN la), and baryon acoustic oscillation (BAO) data. 

We show that there are two CMB shift parameters, R = y^Wfff r(zcMB) (the scaled dis- 
tance to recombination) and l a = ^r(zcMs) It s (zcmb) (the angular scale of the sound horizon 
at recombination), with measured values that are nearly uncorrelated with each other. Allowing 
nonzero cosmic curvature, the three-year WMAP data give R = 1.71 ± 0.03, l a = 302.5 ± 1.2, and 
Q,bh 2 = 0.02173 ± 0.00082, independent of the dark energy model. The corresponding bounds for a 
flat universe are R = 1.70 ± 0.03, l a = 302.2 ± 1.2, and Q b h 2 = 0.022 ± 0.00082. We give the covari- 
ance matrix of (R, l a , 0.bh 2 ) from the three-year WMAP data. We find that (R, l a , Q,bh 2 ) provide 
an efficient and intuitive summary of CMB data as far as dark energy constraints are concerned. 

Assuming the HST prior of H = 72 ± 8 (km/s)Mpc _1 , using 182 SNe la (from the HST/GOODS 
program, the first year Supernova Legacy Survey, and nearby SN la surveys), (R, l a , Qbh 2 ) from 
WMAP three year data, and SDSS measurement of the baryon acoustic oscillation scale, we find 
that dark energy density is consistent with a constant in cosmic time, with marginal deviations 
from a cosmological constant that may reflect current systematic uncertainties or true evolution in 
dark energy. A flat universe is allowed by current data: 0,^ = — 0.006l[5'oi2io 025 f° r assuming that 
the dark energy equation of state wx{z) is constant, and Qk = — 0.002^q'q]J io'oM f° r Wx ( z ) = 
Wo + w a (l — a) (68% and 95% confidence levels). The bounds on cosmic curvature are less stringent 
if dark energy density is allowed to be a free function of cosmic time, and are also dependent on 
the assumption about the early time property of dark energy. We demonstrate this by studying two 
examples. Significant improvement in dark energy and cosmic curvature constraints is expected as 
a result of future dark energy and CMB experiments. 
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I. INTRODUCTION 



The unknown cause for the observed cosmic accelera- 
tion P, Q , dubbed "dark energy" , remains the most com- 
pelling mystery in cosmology today. Dark energy could 
be an unknown energy component 
modification of general relativity 

and [l6| contain reviews of many models. Dark en- 
ergy model-building is a very active research area. For re- 
cent dark ener gy m odels, see for example, [l7l. [l8l. Il9l. I20I 
HI HI HE HI Hi- Current observational data continue 
to be consistent with dark energy being a cosmological 
constant, but the evidence for a cosmological constant is 
not conclusive and more exotic possibilities are still al- 
lowed (see, for example, l2rl [27l . 2p . 

m, m, he in m he \m, m, m 

While the universe is completely consistent with be- 
ing flat under a ACDM hypothesis, it is important to 
note that the observational bounds on dark energy and 
the curvature of the universe are closely related. Cosmic 
Microwave Background (CMB) anisotropy data provide 



the most stringent constraints on cosmic curvature f2fc. 
Assuming that dark energy is a cosmological constant, 
the three-year WMAP data give Q k = -0.15 ± 0.11, 
and this improves dramatically to f2fc = —0.005 ± 0.006 
with the addition of galaxy survey data from the SDSS 
47 1 (2dF data [48( also give a similar improvement) 
6l| . The effect of allowing non-zero curvature on con- 
straining some dark ener gy m odels has been studied by 

he he m, m m, m hem m . 



In this paper, we present a simple and efficient method 
for incorporating constraints from the CMB data into 
an analysis with other cosmological data in constraining 
dark energy without assuming a flat universe. Uisng this 
method, we derive constraints on dark energy and cosmic 
curvature using CMB, type la supernova (SN la) and 
galaxy survey data. 



We describe our method in Sec. II, present our results 
in Sec. Ill, and conclude in Sec. IV. 
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II. METHOD 

The comoving distance from the observer to redshift z 
is given by 

r{z)=cH^ l \n k \- 1 /\mn[\n k \ 1 / 2 T{z% (1) 
r(z) = r^)' E{z)=H{z)/H 

where Qk = —h/H^ with k denoting the curvature con- 
stant, and sinn(x) = sin(x), x, sinh(ir) for Q k < 0, 
fife = 0, and Qk > respectively, and 

e{z) = [n m (i + z f + a-ad(i + zf + n k (i + zf + n x x(z] 

(2) 

with Qx = 1 — fi m — ^rad — £-fc, and the dark energy 
density function X(z) = px(z)/ px(0)- 

CMB data give us the comoving distance to the re- 
combination surface r(zcMs) with zqmb — 1089, and 
the comoving sound horizon at recombination (59l. |60T ] 



o.i 



,1/2 



0.01 



r s (z C MB) 



cH n 



c s dt 
a 

a CMB 



dz 7=7— r 
da 



3(1 + i? b a) a 4 £ 2 (z) 



,(3) 



where a is the cosmic scale factor, clcmb = 1/(1 + 
zcmb), and a A E 2 (z) = Sl m (a + a cq ) + f2 fe a 2 + fix- , -'(^)a 4 , 
with a cq = f2 rad /f2 m = 1/(1 + Zeq)i and z cq = 
2.5 x 10 4 » m /t 2 (r C Mi3 /2.7K)- 4 . The sound speed is 

c s = 1/ -y/ 3(1 + i?6 a), with R b a = 3 i o 6 /(4p 7 ), = 
31500fi b /i 2 (T CM s/2.7K)- 4 . COBE four year data give 
Tcmb = 2.728 K The angular scale of the 

sound horizon at recombination is defined as l a = 

nr{zcMB)/r s {zcMB) 

Note that it is important to use the full expression 
given in Eq.Q in making predictions for l a for dynami- 
cal dark energy models. FigQ] shows how the dark en- 
ergy density X(z) = px{z)/ px{0) compares with the 
matter density p m (z) / p m (0) = (1 + z) 3 for a two pa- 
rameter dark energy model with dark energy equation of 
state wx(z) = Wq + w a (l ~ a) [7(| which corresponds 

to X(z) = a -3(l+«-0+^a) e 3» o (a-l)^ For models with 

wq + w a > 0, the dark energy contribution to the expan- 
sion rate of the universe dominates over that of matter 
at high z. For models that allow significant early dark 
energy (as in the wx(z) — wq + w a (l — a) model), l a 
can be underestimated by 20 — 40% if the dark energy 
contribution to t s (zcmb) is ignored. 1 

We will show that the CMB shift parameters 



R = J fl m H$ r(zcMB), l a = nr(zcMB)/r s (z C MB), 



(4) 
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FIG. 1: Ratio of the dark energy density X(z) = Px { z )/ Px (0) 
and the matter density p m (z)/p m (0) = (1 + z) 3 for dark energy 
models with dark energy equation of state wx (z) = wo + w a (l — a). 



together with ilbh 2 , provide an efficient summary of 
CMB data as far as dark energy constraints go (see 
Sec.IIIA). 

SN la data give the luminosity distance as a function 
of redshift, d L {z) = (1 + z)r(z). We use 182 SNe la 
from the HST/GOODS program [62j and the first year 
SN LS I65| , together with nearby SN la data, as compiled 
by HfTWe do not include the ESSENCE data 
as these are not yet derived using the same method as 
thosed used in [62| • Combining SN la data derived using 
different analysis techniques leads to systematic effects in 
the estimated SN distance moduli [63., [64| . Appendix A 
describes in detail how we use SN la data (flux-averaged 
and marginalized over Hq) in this paper. 

We also use the SDSS baryon acoustic oscillation 
(BAO) scale measurement by adding the following term 
to the x 2 01 a model: 



Xbao 



where A is defined as 



(A - Abao) 



A = 



r 2 (z B Ao) 



czbao 



"I 1/3 



H(ZBAO) 



a 



i H 



1/2 



CZBAO 



(5) 



(6) 



and A bao = 0.469 (ns/0.98) 



-0.35 



OA = 0.017, and 



zbao = 0.35 (independent of a dark energy model) [66j . 
We take the scalar spectral index n s = 0.95 as measured 
by WMAP3 [6l|. 2 



1 The importance of including the dark energy contribution to l a 
is also pointed out by |56H . 



Note that the 16611 constraint on A depends on the scalar spectral 
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For Gaussian distributed measurements, the likelihood 
function L oc e~ x / 2 , with 

X — XCMB + XsNe + XBAO' V< ) 

where Xcmb * s gi ven m Eq. (TiT))) in Sec.IIIA, XsNe 1S gi ven 
in Eq. (|A3p in Appendix A, and Xbao 1S g iven m Eq.([5]). 

We derive constraints on the dark energy density func- 
tion X(z) = px{z) / pxifi) as a free function at z < z cut , 
with its value at redshifts z% = z cu t{i/n) (i=l, 2, n), 
X{zi), treated as n independent parameters estimated 
from data. We use n = 3 and z cut = 1.4 in this paper. 
We use cubic spline interpolation to obtain values ofX(z) 
at other values of z at z < z cut |26j . The number of cur- 
rently published SNe la is very few beyond z cut = 1.4. 
For z > z cu t, we assume X(z) to be matched on to either 
a powerlaw [26| : 

X(z)=X(z cut )(-^-Y , (8) 

or an exponential function: 

X(z)^X(z cut )e a ^ zcut \ (9) 

We impose a prior of a > —3 as a is not bounded from 
below. Our approach effectively decouples late time dark 
energy (which is responsible for the observed recent cos- 
mic acceleration and is probed directly by SN la data) 
and early time dark energy (which is poorly constrained) 
by parametrizing the latter with an additional parameter 
estimated from data. 

For comparison with the results of others, we also de- 
rive constraints for models with dark energy equation of 
state wx(z) = wq + w a (l — a). This parametrization has 
the advantage of not requiring a cutoff to obtain a finite 
dark energy equation of state at high z (which is not 
true for the wx(z) — Wq + w'z parametrization), but it 
does allow significant early dark energy (which can cause 
problems for Big Bang Nucleosynthesis [87J 3 and cos- 
mic structure formation [8Sj). unless a cutoff is imposed. 
This dilemma illustrates the limited usefulness of simple 
parametrizations of dark energy. 

For all the dark energy constraints from combining the 
different data sets presented in this paper, we marginal- 
ize the SN la data over Hq in flux-averaging statistics 



index ng. Since the error on nj from WMAP data does not 
increase the effective error on A, and the correlation of rig with 
R, and l a is weak, we have ignored the very weak correlation of A 
with _R and l a in our likelihood analysis. We have derived R and 
I a from WMAP data marginalized over all relevant parameters. 
3 The current BBN constraints, S = 0.942 ± 0.030 (N v = 
2.30^0 34) rule out the standard model of particle physics (S = 1, 
N v = 3) at la [87ll , Given the uncertainties involved in deriving 
the BBN constraints, we relax the standard deviation of S by a 
factor of two, so that the standard model of particle physics is 
allowed at ltr. Wc find that the resultant BBN constraints do 
not have measurable effect on our dark energy constraints. 



(described in the next subsection) , and impose a prior of 
H = 72 ± 8 (km/s)Mpc- 1 from the HST Cepheid vari- 
able star observations fn\ . 

We run a Monte Carlo Markov Chain (MCMC) based 
on the MCMC engine of ^ to obtain O(10 6 ) samples 
for each set of results presented in this paper. For the 
full CMB analysis we used the WMAP three year tem- 
perature and polarization 4 power spectra [6l| with ver- 
sion 2 of their likelihood code [68| together with theo- 
retical power spectra generated by CAMB (with pertur- 
bations in dark energy) [69j ; the parameters used are 
(fife, fib/1 2 , £! m /i 2 , h, A s , r, n s , p DE ). For the com- 
bined data analysis using CMB shift parameters, the pa- 
rameters used are (fife, fi m , h, fif,/i 2 , Pde)- The dark 
energy parameter set Pde = w for a constant wx(z), 
Pde = (wo,w a ) for w x {z) = w + w a (l - a), and 
Pde = (X (zi) , X (z 2 ) , X (z 3 ) , a) for the general case. We 
assumed flat priors for all the parameters, and allowed 
ranges of the parameters wide enough such that further 
increasing the allowed ranges has no impact on the re- 
sults (with the exception of constraining w and (u>o, w a ) 
using CMB data only where we have to impose fixed al- 
lowed ranges for w and (wo,w a ) since these are not well 
constrained). The chains typically have worst e- values 
(the variance(mean)/mean(variance) of 1/2 chains) much 
smaller than 0.01, indicating convergence. The chains 
are subsequently appropriately thinned to ensure inde- 
pendent samples. 



III. RESULTS 

A. A Simple and Efficient Method for 
Incorporating CMB data 

1. A roadmap of our method 

We propose a simple and efficient method for dark 
energy data analysis, with \ 2 — ~ 2 hi£ = Xcmb + 
XsNe + XbaOi where Xcmb is g iven b Y constraints on 
{R,l a ,Vl b h 2 ) (sec Eq.[TU] in Sec.IIIA), x"sNe ^ s given by 
SN la data flux-averaged and marginalized over Hn (see 
Eq. |A3j in Appendix A) , and Xbao ^ s gi ven by [6(| (see 
Eq. [S])- In our method, CMB data are incorporated by 
using constraints on (_R, l ai fi^/i ), instead of using the 
full CMB power spectra. In Sec lIII A2l below, we will 
show that (R, l a , fi^/i 2 ) provide an efficient and intuitive 
summary of CMB data as far as dark energy constraints 
are concerned. 



4 The main contribution of CMB polarization data is the determi- 
nation of the reionization optical depth. 
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2. Justification of our method 

We have performed MCMC calculations using only the 
full CMB temperature and polarization angular power 
spectra from WMAP three year observations, without 
assuming spatial flatness, and without imposing any pri- 
ors on H . These calculations are quite time consuming. 
We have used these to derived the results in Fig[2] and 
Tables I II. 

FigO shows that allowing nonzero cosmic curvature, 
the three-year WMAP data give measurements of (R, l a , 
Vlbh 2 ) that are independent of the dark energy model. 5 
The measurements of (R, l a , f^ft 2 ) differ slightly in 
a flat universe because of the correlation of curvature 
with other cosmological parameters when spatial flatness 
is not assumed. Table I gives the parameters for the 
Gaussian fits to the probability distribution functions of 
(RJ a ,tt b h 2 ,n m h 2 ,r s (z CMB ),r(z CMB )) from the three- 
year WMAP data. 6 These fits are independent of the 
dark energy model assumed. The constraints on (fl m h 2 , 
t s {zcmb), r(zcMB)) are also independent of the assump- 
tion about cosmic curvature. 

Table II gives the normalized covariance matrices for 
(R,la,ftbh 2 ,£l m h 2 ,r s (zcMB),r(zcMB)) from the three- 
year WMAP data for a ACDM model for models with 
and without curvature. These are appropriate to use with 
Table I; models with non-constant dark energy density 
give slightly smaller correlations between the parameters. 
Note that we have included (fl m h 2 , t s {zcmb), ^{zcmb)) 
in Tables I-II to show that although these three parame- 
ters are well constrained by CMB data, they are strongly 
correlated with each other 7 , in contrast to the parameters 
we have chosen, (R, l a , f2f,/i 2 ). 

We find that there are two CMB shift parameters, R 
and l a (with measured values that are nearly uncorre- 
cted, see Table II), that are optimal for use in constrain- 
ing dark energy models. 8 FigJ3] shows that both R and 



5 R. and l a are shifted slightly if the running of ng or/and a nonzero 
tensor to scalar ratio are considered, and shifted more notably 
if a nonzero neutrino mass is considered l8ll|. C urrent CMB data 
do not require these additional parameters [6l||. 

6 Note that CMB data do not constrain Ho in models with nonzero 
curvature due to parameter degeneracies. For example, the di- 
mensionless Hubble constant h = 0.50 ±0.14 for a ACDM model 
with Hj, 7^ 0. It is the absolute scales of t s (zcmb) and t(zcmb) 
that are well determined by the CMB data. 

7 This high degree of correlation arises from how these three pa- 
rameters are measured. The sound horizon at recombination 
r s {zc M b) is derived primarily using the measurements of C(,?t 2 
and Q m h 2 [60H . hence is strongly correlated with fl m h 2 . The 
distance to the recombination surface rjzcMfl) is derived using 
t s{zcmb) an d the angular scale of the sound horizon l a [60|. 
hence is strongly correlated with t s (zcmb) l60ll 

8 R has been known as the CMB shift parameter in the past 
[Hi IM IzH ■ El showed that in an open universe with a 
cosmological constant, there is a degeneracy along the SR = 
lines, i.e., models with different values of fl m , ^a, and h that 
give the same value of R are not distinguishable except at very 
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FIG. 2: The scaled distance to recombination R, the angular scale 
of the sound horizon at recombination l a , and the baryon density 
Ufjh 2 from the three- year WMAP data. 
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TABLE I: The parameters for the Gaussian fits to the probability distribution functions of 
(R, l a , tibh 2 , Sl m /i 2 , r s (zcMB), r(zcMB)) from the three-year WMAP data, independent of the dark energy model assumed. 



Parameter 


mean 


rms variance 


^ m h 


0.1284 


0.0086 


r s (z G MB) /Mpc 148.55 


2.60 


r(zcjffl)/Mpc 


14305 


285 




R 


1.71 


0.03 


la 


302.5 


1.2 


Q b h 2 


0.02173 


0.00082 




R 


1.70 


0.03 


la 


302.2 


1.2 


Q b h 2 


0.022 


0.00082 



TABLE II: Normalized covariance matrices for (R, l a , Q b h 2 , ft m h 2 , r s (zcMB), t{zcmb)) fromthe WMAP three year data. 





R 


la 


n b h 2 


^ m h 


r s {zcMB) 


r(zcMB) 


n k /o 


R 


0.1000E+01 


-0.1237E+00 


0.6627E-01 


0.9332E+00 


-0.8805E+00 


-0.8023E+00 


la 


-0.1237E+00 


0.1000E+01 


-0.6722E+00 


-0.4458E+00 


0.5214E+00 


0.6569E+00 


n b h 2 


0.6627E-01 


-0.6722E+00 


0.1000E+01 


0.3731E+00 


-0.5047E+00 


-0.5778E+00 




0.9332E+00 


-0.4458E+00 


0.3731E+00 


0.1000E+01 


-0.9882E+00 


-0.9605E+00 


Ts(zcMb) 


-0.8805E+00 


0.5214E+00 


-0.5047E+00 


-0.9882E+00 


0.1000E+01 


0.9859E+00 


t{zcub) 


-0.8023E+00 


0.6569E+00 


-0.5778E+00 


-0.9605E+00 


0.9859E+00 


0.1000E+01 


n k = 


R 


0.1000E+01 


-0.904 7E-01 


-0.1970E-01 


0.9397E+00 


-0.8864E+00 


-0.8096E+00 


la 


-0.9047E-01 


0.1000E+01 


-0.6283E+00 


-0.3992E+00 


0.4763E+00 


0.6185E+00 


n b h 2 


-0.1970E-01 


-0.6283E+00 


0.1000E+01 


0.2741E+00 


-0.4173E+00 


-0.4942E+00 




0.9397E+00 


-0.3992E+00 


0.2741E+00 


0.1000E+01 


-0.9876E+00 


-0.9594E+00 


r s (zcMB) 


-0.8864E+00 


0.4763E+00 


-0.4173E+00 


-0.9876E+00 


0.1000E+01 


0.9855E+00 


t(zcmb) 


-0.8096E+00 


0.6185E+00 


-0.4942E+00 


-0.9594E+00 


0.9855E+00 


0.1000E+01 



l a must be used to describe the complex degeneracies 
amongst the cosmological parameters that determine the 
CMB angular power spectrum. 

Fig[3] illustrates the relationship of R and l a in deter- 
mining the CMB angular power spectra for simple models 
that give the same R or l a values. Fig[3{a) shows that 
models that correspond to the same value of R but differ- 
ent values of l a give rise to very different CMB angular 
power spectra because l a determines the acoustic peak 
structure. Fig[3]Jb) shows that models that correspond 
to the same value of l a but different values of R have 



the same acoustic peak structure in their CMB angular 
power spectra, but the overall amplitude of the acoustic 
peaks is different in each model because of the difference 
in R. 9 

Now we illustrate how using both R and l a helps con- 
strain models with a constant dark energy equation of 
state, and zero or small curvature (the class of models 
shown in Fig[3]). FigfJ] shows the expected R, l a and 
Ts(zcmb) as functions of fl m for five models. For refer- 
ence, the values for h and flf,h 2 have been chosen such 
that the cosmological constant model satisfies both the 



low multiples (where cosmic variance dominates), see Fig.l of 
their paper. 



R is proportional to Cl m h 2 , which determines the overall height 
of the acoustic peaks. 
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FIG. 3: CMB angular power spectra for dark energy models that 
give the same values of R or l a . 



FIG. 4: The expected R and l a as functions of curvature for five 
simple dark energy models (with the same line types as in Fig[3j . 



R and l a constraints from WMAP three year data at the 
same value of f2 m (as in Fig|3]). Note that for the other 
four models, the R and l a constraints cannot be satisfied 
at the same Q m value. This is because R and t s (zcmb) 
have different dependences on f2 m . Models that give the 
wrong R and r s (zcMB) values can give the right value for 
l a because l a oc R/r s (zcMB)- Using both R and l a con- 
straints thus helps tighten the constraint on f2 m , which 
leads to tightened constraints on w or f^. 

When more complicated dark energy models and 
nonzero cosmic curvature are considered, there is a de- 
generacy between dark energy density function X(z) and 
curvature. The R or l a constraints from CMB can al- 
ways be satisfied with a suitable choice of curvature, 
but satisfying the R and the l a constraints usually re- 
quire different values for curvature. Thus using both R 
and l a constraints from CMB helps break the degeneracy 
between dark energy parameters and curvature. Fig[5] 
demonstrates this by showing the expected R, l a , and 
t 8 {zcmb) as functions of curvature for the dark energy 
models from FigfT] (with the same line types). For refer- 
ence, the values for Vt m and h have been chosen such that 
the cosmological constant model satisfies both the R and 
l a constraints from WMAP three year data. Clearly, the 
R constraint rules out closed models with large curva- 



ture, while the l a constraint rules out open models with 
large curvature. The vertical dotted lines indicate the 
1 a range of fl^ from R, l a , and f2&/i 2 constraints from 
WMAP three year data, combined with the data of 182 
SNe la, and the SDSS BAO measurement. 

Note that the baryon density ilbh 2 should be included 
as an estimated parameter in the data analysis. This is 
because the value of fl^h 2 is required in making a predic- 
tion for l a in a given dark energy model (see Eq. [3]), and 
it is correlated with l a (see Table II). 

To summarize, we recommend that the covariance ma- 
trix of (R, l a , Q,iyh?) given in Tables. I-II be used in the 
data analysis. To implement this, simply add the follow- 
ing term to the \ 2 of a given model with pi = R, y% = l a , 
and P3 = fli,h 2 : 



Xcmb = A Pi l Cov 1 (puPj)] A Pj> 



Pi -Pi 



data 



(10) 



where pf ata are the mean values given in Table I. The 
covariance matrix Cov(pi,pj) is obtained by multiply- 
ing the normalized covariance matrix in Table II with 
[Var(pi) Varfjjj)] 1 / 2 , with the rms variance [Var(pi)] 1 / 2 
given in Table I. Note that our constraints on (R, l a , 
flbh 2 ) have been marginalized over all other parameters 
including the dark energy parameters. 

As a test for the effectiveness of our simple method for 
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B. Constraints on dark energy 
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FIG. 5: The expected R and l a as functions of curvature for the 
dark energy models from Fig[T](with the same line types). Both R 
and l a are needed to constrain cosmic curvature. 



incorporating CMB data, we derived the constraints on 
wx(z) = w (constant) and u>x{z) = WQ + w a (l — a) using 
(i?, l a , flbh 2 ), and compared with the results from using 
the full CMB code CAMB. For both sets of calculations, 
we assumed the same flat priors of— 2 < w < 0, — 2 < 
wq < 0, and —6 < w a < 3, since w and (wo,w a ) are not 
well constrained by using CMB data alone. The pdf's of 
w and (wo,w a ) span the entire allowed ranges, and have 
similar shapes in the two methods. We did not assume 
any priors on Hq since we want to study CMB data only. 
For wx(z) = w (constant), using (R, l a , Q^h 2 ) gives w = 
-0.96 ±0.57, while the full CMB code CAMB gives w = 
-0.97 ±0.53. Using (i?, l a , Vl b h 2 ) w a = -1.0 ±0.6 and 
w a = -2.2 ± 2.1, while the full CMB code CAMB gives 
wo = — 0.9±0.6 and w a = — 2.4±2.0. These comparisons 
indicate that our simple method of incorporating CMB 
data by using Eq. (|10[) is indeed efficient and appropriate 
as far as dark energy constraints are concerned. Since 
CMB data alone do not place tight constraints on dark 
energy, it is not appropriate to do the comparison of our 
method with the full CMB code for dark energy models 
with more parameters. 



Because of our ignorance of the nature of dark energy, 
it is important to make model-independent constraints 
by measuring the dark energy density px(z) as a free 
function. Measuring px(z) has advantages over measur- 
ing dark energy equation of state wx(z) as a free func- 
tion; px{z) is more closely related to observables, hence 
is more tightly constrained for the same number of red- 
shift bins used [76|, [77], [zH . More importantly, measuring 
wx(z) implicitly assumes that px(z) does not change 
sign in cosmic time (as px{z) is given by the exponen- 
tial of an integral over 1 + wx(z)); this precludes whole 
classes of dark energy models in which px{z) becomes 
negative in the future ("Big Crunch" models, see [89[ for 
an example) [26[ . 

We have reconstructed the dark energy density func- 
tion X(z) = px{z)/ px(0) by measuring its value at 
Zi = z cut (i/3) (i=l, 2, 3) at z < z cut , and parametrized 
it by either a powerlaw (X(z) oc (1 + z) a ) or an exponen- 
tial function (X(z) oc e az ) at z > z cut (see Eqs.dH])-®). 
We have chosen z cut — 1.4 as few SNe la have been ob- 
served beyond this redshift. We find that current data 
allow a > for X{z) oc (1 + z) a at z > z cut , and re- 
quire a < for X{z) oc e az at z > z cut . This means that 
assuming powerlaw dark energy at early times allows sig- 
nificant amount of dark energy at 2 > 1, while assuming 
exponential dark energy at early times is equivalent to 
postulating dark energy that disappears at z 1. The 
latter is more physically sensible since dark energy is in- 
troduced to explain late time cosmic acceleration. Intro- 
ducing dark energy that is important at early times could 
cause problems with Big Bang Nucleosynthesis [87j and 
formation of cosmic large scale structure |88[ . 

Fig[S] shows the reconstructed dark energy density 
function X(z) using (R, l a , f2(,/i 2 ) from the three-year 
WMAP data, together with 182 SNe la and SDSS BAO 
measurement. The apparent shrinking of the error con- 
tours at z > z cut is due to the use of one parameter 
to describe X(z) at z > z cut . Future theoretical work 
and better data will allow better-motivated description 
of dark energy at early times. 10 FigJT] shows the cor- 
responding constraints on the cosmic expansion history 
H(z). 

For a flat universe, the dark energy constraints at 
z < 1 are nearly independent of the early time assump- 
tion about dark energy, while the dark energy constraint 
at z ~ z cut is more stringent if X(z) oc (1 + z) a at 
z > z cut . This is as expected. Because of parameter 
correlations, stronger assumption about early time dark 
energy (the powerlaw form) leads to more stringent dark 
energy constraint at late times around z ~ z cut . 

Without assuming a flat universe, in the X(z) oc 
(1 + z) a at z > z cut case, there is a strong degeneracy 



See for example, 19011 . which assumed a flat universe. 
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FIG. 6: Constraints on dark energy density using R, l a , and Qj,/i 2 
from the three-year WMAP data, together with 182 SNe la and 
SDSS BAO measurement. The shaded areas indicate the 68% con- 
fidence regions, while the outside contours bound the 95% confi- 
dence regions. 



FIG. 7: Constraints on the expansion history of the universe H(z) 
that corresponds to Fig[6] using Ft, l a , and flf,h 2 from the three- 
year WMAP data, together with 182 SNe la and SDSS BAO mea- 
surement. The error bars indicate the 68% confidence intervals. 



between curvature and the powerlaw index a. This is 
as expected since the curvature contribution to the total 
matter-energy density is also a powerlaw, (1 + z) 2 . X(z) 
is not well constrained in this case, and is not shown in 
FigJS] When X(z) oc e az is assumed at z > z cut , there is 
no degeneracy between the exponential index a and cur- 
vature. X(z) is well constrained in this case (see Fig|6]). 

For comparison with the work by others, Fig[8] shows 
the constraints on (wo,w a ) for models with dark energy 
equation of state u>x(z) = Wq + w a (l — a), using R, l a , 
and £!{,/i 2 from the three-year WMAP data, together with 
182 SNe la and SDSS BAO measurement. These are 
consistent with the results of [56|, [5?| • Note that using 
wx (z) = wo + w a ( 1 — a) implies extrapolation of dark en- 
ergy to early times, which leads to artificially strong con- 
straints (compared to model-independent constraints) on 
dark energy at both early and late times. This was noted 
by [fill as well. 

Comparing Fig[6][8]with Figs. 3-5 of [75[ (for the case of 
assuming X{z) oc (l + z) a at z > z cu t), it is clear that the 
constraints on dark energy have significantly tightened if 
a flat universe is assumed. 



C. Cosmic curvature and dark energy constraints 

Fig[9] shows the probability distribution function of 
cosmic curvature for different assumptions about dark en- 
ergy: the model- independent dark energy density px(z) 
reconstructed in the last subsection, the two param- 
eter dark energy model wx(z) — wq + w a (l — a), 
and a constant dark energy equation of state. A flat 
universe is allowed at the 68% confidence level in all 
the cases when curvature is well constrained, fifc = 
— 0.006+°; 012^0 025 f° r assuming that w x {z) is constant, 
and n k = -0.0021^1°;°^ for w x (z) = w +w a (l- a) 
(68% and 95% confidence levels). Assuming a constant 
dark energy equation of state gives the most stringent 
constraints on cosmic curvature. The bounds on cosmic 
curvature are less stringent if dark energy density is al- 
lowed to be a free function of redshift, and are dependent 
on the assumption about the early time property of dark 
energy. If dark energy is assumed to be an exponential 
function at z > z cut (z cut = 1.4), it is well constrained 
by current observational data (see FigO and negligible 
at early times. In this case, curvature is well constrained 
as well. If dark energy is assumed to be a powerlaw at 
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FIG. 8: Constraints on (wo,w a ) using R, l a , and Q[,/i 2 from the 
three-year WMAP data, together with 182 SNe la and SDSS BAO 
measurement. The 68% and 95% confidence contours are shown. 



early times, its powerlaw index is strongly degenerate 
with curvature, and neither is well constrained. 



IV. SUMMARY AND DISCUSSION 

We have presented a simple and effective method for 
incorporating constraints from CMB data into an analy- 
sis of other cosmological data (for example, SNe la and 
galaxy survey data), when constraining dark energy with- 
out assuming a flat universe. 

We find that three-year WMAP data give constraints 
on (R,l a ,n b h 2 ,£l m h 2 ,r s (zcMB),r{z C MB)) that are in- 
dependent of the assumption about dark energy (see Ta- 
ble I). The constraints on (fi m /i 2 , t s (z C mb), r{zcMB))) 
are also independent of the assumption about cosmic cur- 
vature, but they are strongly correlated with each other 
and are not suitable for use in constraining dark energy 
(see Table II). 

We show th at there are two CMB shift parameters, 
R = yjQ m H'g t(zcmb) (the scaled distance to recombi- 
nation) and l a = itrizcMB) I t s (zcmb) (the angular scale 
of the sound horizon at recombination); these retain the 
sensitivity to dark energy and curvature of t(zcmb) and 



1.2 r i i i i i 

: WMAP 3 yr + 182 SNe la + SDSS BAO 




FIG. 9: The probability distribution function of cosmic curva- 
ture for different assumptions about dark energy: the model- 
independent dark energy density px( z ) reconstructed in the last 
subsection, the two parameter dark energy model wx{z) = wq + 
ui a (l — a), and a constant dark energy equation of state. 



t s {zcmb), and have measured values that are nearly un- 
corrected with each other (see Table II) . We give the co- 
variance matrix of (R, l a , Slbh 2 ) from the WMAP three 
year data (see Tables I and II). 

We demonstrate that (R, l a , £lbh 2 ) provide an efficient 
summary of CMB data as far as dark energy constraints 
are concerned, and an intuitive way of understanding 
what the CMB does in terms of parameter constraints 
(see FigsHS]). 

While completing our paper (based on detailed calcula- 
tions that have taken several months), we became aware 
of Ref. 81]. They also found that using both R and l a 
tightens dark energy constraints. However, their paper 
assumed a fiat universe, and used an approximation for 
l a that ignores both curvature and dark energy contribu- 
tions. We use the exact expression for l a and derived the 
covariance matrix for (i?, l a , fl^h 2 ) which are based on 
the MCMC chains from our full CMB power spectrum 
calculations without assuming spatial flatness. 

We have used (R, l a , Q b h 2 ) from WMAP three year 
data, together with 182 SNe la (from the HST/GOODS 
program, the first year Supernova Legacy Survey, and 
nearby SN la surveys), and SDSS measurement of the 
baryon acoustic oscillation scale in deriving constraints 
on dark energy. Assuming the HST prior of Hq = 
72 ± 8 (km/s)Mpc _1 [7l|, we find that current observa- 
tional data provide significantly tightened constraints on 
dark energy models in a flat universe, and less stringent 
constraints on dark energy without assuming spatial flat- 
ness (see Figs|B][5]). Dark energy density is consistent 
with a constant in cosmic time, with marginal deviations 
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from a cosmological constant that may reflect current 
systematic uncertainties 11 or true evolution in dark en- 
er gy ( see Fi gs EH ). Our findings are consistent with that 
of [62[ and [431 ]. 

A flat universe is allowed by current data at the 68% 
confidence level. As expected, the bounds on cosmic cur- 
vature are less stringent if dark energy density is allowed 
to be a free function of cosmic time, and are also de- 
pendent on assumption about dark energy properties at 
early times (see FigEJ) - The behavior of dark energy at 
late times (where it causes cosmic acceleration and is di- 
rectly probed by SN la data) and at early times (where it 
is poorly constrained) should be separated in parameter 
estimation in order to place robust constraints on dark 
ener gy a nd cosmic curvature (see Sec.IIIB and C). 

11 Ref l37H studied the statistical consistency of subsets of SNe la 
that comprise the 182 SNe la. 



Future dark energy experiments from both ground 
and space (8 2. pot [8J, [85]], together with CMB data 
from Planck [86|, will dramatically improve our ability 
to probe dark energy, and eventually shed light on the 
nature of dark energy. 
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APPENDIX A: MARGIN ALIZATION OVER Ho 
IN SN IA FLUX STATISTICS 

Because of calibration uncertainties, SN la data need 
to be marginalized over Hq if SN la data are combined 
with data that are sensitive to the value of Hq. This is 
the case here (see the next section). We use the angu- 
lar scale of the sound horizon at recombination l a which 
depends on £l m h 2 , while the dimcnsionless Hubble pa- 
rameter E(z) = H(z)/Hq (which appears in the deriva- 
tion of all distance-redshift relations) depends on O m . 
Hence a dependence on Hq is implied. We marginal- 
ize the SN la data over Hq while imposing a prior of 
Hq = 72 ± 8 (km/sftlpc- 1 from HST Cepheid varibale 
star observations [7ll |. 

The marginalization of SN la data over Hq was derived 
in [76| for the usual magnitude statistics (assuming that 
the intrinsic dispersion in SN la peak brightness is Gaus- 
sian in magnitudes). Here we present the formalism for 
marginalizing SN la data over Hq in the flux-averaging of 
SN la data using flux statistics (see Eq. |Al| ). The pub- 
lic software for implementing SN la flux averaging with 
marginalization over Hq (compatible with cosmomc) is 
available at http:/ /www. nh n.ou .e du/~w ang/SNc ode/| 

Flux-averaging of SN la data [6J| is needed to minimize 
the systematic effect of weak lensing of SNe la [7!| . [lH 
presented a consistent framework for flux-averaging SN 
la data using flux statistics. Normally distributed mea- 
surement errors are required if the x 2 parameter estimate 
is to be a maximum likelihood estimator 80]. Hence, if 
the intrinsic dispersion in SN la peak brightness is Gaus- 
sian in flux, we have 



XN da 



E 



[F( Zl )-FV(z t \s)Y 



(Al) 



Since the peak brightness of SNe la have been given in 
magnitudes with symmetric error bars, m pea k ± Cm, we 
obtain equivalent errors in flux: 



After flux-averaging, we have 



^ \F(zj)- FP(zj js) 

i FA 



(A2) 



where f(Z«|si) = (d L (z l |s)/Mpc) \ 

The predicted SN la flux FP( Zi \s) = 
[di(v2i|s)/Mpc] 2 oc h 2 . Assuming that the dimcn- 
sionless Hubble parameter h is uniformly distributed in 
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the range [0,1], it is straightforward to integrate over h , , _ , 2 .^v^ ^( z i) (aa\ 

in the probability distribution function to obtain yo\ ) — \ ) ^a— ' 

i F,i 



P (s\0 <h<l) = e"*^ = Jo*™ (A3) 

J dxe~9°W 



where 



c/(a;) = 2^ 5-2 — . where F*{Zi\s) = FP(zi\s, h = 1). 



